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ABSTRACT 



Protostellar outflows crisscross the regions of star cluster formation, stirring turbulence and 
altering the evolution of the forming cluster. We model the stirring of turbulent motions by 
protostellar outflows, building on an observation that the scaling law of supersonic turbulence 
implies a momentum cascade analogous to the energy cascade in Kolmogorov turbulence. We 
then generalize this model to account for a diversity of outflow strengths, and for outflow col- 
limation, both of which enhance turbulence. For a single value of its coupling coefficient the 
model is consistent with turbulence simulations by Li & Nakamura and, plausibly, with obser- 
vations of the NGC 1333 cluster-forming region. Outflow-driven turbulence is strong enough 
to stall collapse in cluster-forming regions for several crossing times, relieving the mismatch 
between star formation and turbulent decay rates. The predicted line-width-size scaling implies 
radial density indices between -1 and -2 for regions supported by outflow-driven turbulence, 
with a tendency for steeper profiles in regions that are more massive or have higher column 
densities. 



1. INTRODUCTION 

Star formation in the Milky Way is known to occur 
predominantly within star clusters. Clusters with dif- 
ferent multiplicites are born in a remarkably flat 
distribution, dAf c \/dln A* oc . This is flat in the 
sense that the birth cohort of a given star is equally 
likely to be in any decade of JV* within the allowed 
range, from a lower limit of A* ~ 50—100 to an upper 
limit that depends on the area sampled - a few thou- 



dense) context of the Sun's birth cluster. 

There is growing evidence that star cluster forma- 
tion is a slow process, in the sense that the star forma- 
tion rate, though quite uncertain, is about an order 
of magnitude slower tha n the mean free-fall ra te of 
the parent gas "clump". Palla h Stahler ( 19991 ) and 
Huff fc Stabled ([2006) argue from the Orion Nebula 



sand , based on surveys out to 2 kpc (ILada Sz L ada 
2003J), or about half a million (IMcKee &: Williams 
19971 ). based on HII region surveys of the entire in- 
ner Galaxy. The number, density, and lifetime of 
the birth cohort affect the likelihood of exposure to 
ionizing radiation and supernovae from the massive 
cluster members, as well as the likelihood of close stel- 
lar encounters, all of which affect the formation and 
evolution of planetary systems. 

The Solar System, for instance, is thought to 
deriv e from a cluster with a few thousand mem- 
bers ( Adams Sz Laughlinl l200ll ) on the basis that a 
nearby supernova supplied short-lived radionuclides, 
and that passing stars did not disturb t he outer plan- 
ets. The orbit of the planetoid Sedna (IBrown et al 



Cluster's Hertzsprung-Russell diagram that star for- 
mation began there at least 10 Myr ago, as has been 
dramatically confirmed by the de tection of lithium 
depletion in some of its members ( Palla et al. 20051 ). 
The common origin of runaway stars AE Aurigae 
and n Columbae in the clu ster's vicinity 2.6 Myr 
ago (IBlaauw k, Morganlll954r) supports this assertion 
(jTan et alJl200d ). I Tan et al.l also point to the fact 
that cluster mass profiles are quit e smooth, requir- 
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2005) appears to require an encounter with another 
star well within a thousand AU, an event which is 
nearly impossible except in the dense (but not too 



i ng tim e for subclusters to merge. 
(2006) compute the star formation rate parameter 
SFRff , which is normalized to the free-fall rate, from 
molecular and infrared observations of dense, cluster- 
forming clumps, finding SFRff ~ 10~ 17 (with signifi- 
cant uncertainty) for hydrogen densities ranging from 
riff = 10 2 to 10 53 cm~ 3 . As they point out, this is 
consistent both with theories in which star forma- 
tion is regul ated by ambip olar diffusion within mag- 
netized gas ( McKed Il989i ). and with those in which 
turbulence suppresses the rate of localized collapse 
(jKrumholz & McKed 120051 ) . Regardless of its cause, 
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slow cluster formation is in sharp contrast with the 
rapid decay of superso nic turbulence over at most 
a cou ple free-fall times ( Stone et al.l 19981 ; Mac Low 
1 19991 ). Therefore, forming star clusters must tap 
some source of energy to sustain the observed tur- 
bulence - and the nature of this feedback may leave 
its imprint on the clusters' properties. 

This paper addresses the driving of turbulence by 
protostellar jets and winds, which are thought to 
emerge from stars of all masses as they form. Since 
outflows driven by these jets and winds are ubiqui- 
tous, outflow-driven turbulence is a pervasive and in- 
evitable consequence of stellar cluster formation. 

A second, more sporadic form of feedback comes 
from the most massive cluster members, whose 
ionizing radiation is a strong driver of turbulence 
within giant molecular clouds (GMCs; iMatznei 2002; 
iKrumholz et al. I l200d ). However, cluster-forming 
clumps are much denser than the GMCs surrounding 
them (~ 10 4-5 rather than ~ 10 2 H atoms cm -3 ), 
and high densities damp the dynamical effects of 
photoionization. We shall ther efore concentrate on 
outflow-driven turbulence. (See ITan Sz McKee 2000 
for preliminary cluster-formation models which in- 
cl ude photoioniza ti on.) 

INorman h Silk! ( 1980l ) proposed that protostel- 
lar outflows sustain the turbulence in GMCs, be- 
fore it was known that star formation is re- 
stricted to the densest and m ost massive clumps 
within them. McKee ( 19891 ) adopted protostel- 



lar outflows as the driving agents in his equi- 
librium theory of GMCs. Conce ntrating o n the 
scale of a cluster- form in g clum p, iMatzneil (1999) 
and Matzner Sz McKee ( 1999bl ) incorporated out- 
flows into dynamical models for star cluster forma- 
tion. Recent s i mula tions of cluster formation by 
Li Nakamural (|2006l ) provide an important demon- 



stration that protostellar outflows are dynamically 
important in this process. 

The theory of outflow-driven turbulence presented 
here improves upon these previous studies by ac- 
counting for two of the outflows' key properties: a 
strong collimation toward an outflow axis, and a di- 
versity in strength associated with the range of pos- 
sible stellar masses. These effects conspire to en- 
hance the outflows' effect on the dynamics of cluster- 
forming clumps, relative to a model in which outflows 
are spherical and all have the same strength. 

In §[2] we establish the basic dimensional scales and 
dimensionless ratios that define outflow-driven turbu- 
lence. In £§] we propose models for the energy spec- 
trum and line-width-size relation of outflow-driven 
turbulence, drawing on an insight about the dynami- 
cal nature of supersonic turbulence in strongly radia- 
tive gas. We extend this model to a diversity of out- 



flow intensities in § 14. H and evaluate it for the stellar 
initial mass function (IMF) and outflow collimation 
in §4.1.11 and 14.1.21 respectively. The loss of momen- 
tum from a finite region due to collimated outflows 
is considered in § 14.21 We compare to available sim- 
ulations in §[5] and to observations of NGC 1333 in 
§E1 before drawing conclusions about the dynamics 
of stellar cluster formation in §0 

2. DIMENSIONAL ANALYSIS 

To begin as generally as possible, let us first con- 
sider a pressureless medium of mean density po stirred 
by spherical explosions whose characteristic momen- 
tum is I, that occur at a rate per unit volume S. We 
choose to describe these by their momentum rather 
than their energy, because we expect the shocked ma- 
terial to be highly radiative. These three parameters 
with independent dimensions define scales of mass, 
length, and time: 
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One quantity of interest is the turbulent line-width, 
i.e., the characteristic turbulent velocity, measured 
on scale £: 

(2) 



t 



4/7 

Po 



10 39,6 gcms 1 
10" 19 - 6 gem" 3 



The dynamical state of the gas will be highly tur- 
bulent and clumpy, but in the absence of other di- 
mensional parameters any quantity must obey these 
scalings. There is one explosion, on average, per £ 3 t. 
Intermittent explosions are the dominant feature of 
flow on scales smaller than I, whereas flow on larger 
scales reflects many overlapping explosions. 

To estimate characteristic scales for the cluster- 
forming regions discussed later, consider I = 
S = 10~ 67 - 2 cm" 3 s" 1 , and p = 
values typical of cluster-forming re- 
gions like NGC 1333 (§ EJ). Then, £ = 0.38 pc, 
t = 0.34 Myr, m = 19 M , and £/t = 1.1 kms" 1 . 
Outflow-driven turbulence (1) extends well beyond 
the regions of localized collapse associated with in- 
dividual forming stars, and (2) involves motions sig- 
nificantly faster than the thermal sound speed (c s ~ 
0.2 kms" 1 ). In § !4.1l we shall find that these qualities, 
which distinguish outflows as likely agents of support, 
are only enhanced by effects like outflow collimation. 

2.1. Gravity, stellar mass, and finite-size regions 

We now introduce gravity, and restrict our atten- 
tion to a region of finite size R - hence of mass 
M = AttB?pq/J> if we take po to be the mean den- 
sity. Gravity and finite radius come together nat- 
urally, because star formation is known to occur in 
self- gravitating "clumps" of finite mass and radius. 
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Let us also introduce a stellar mass scale M*, and 
suppose that the explosion momentum is related to 
M* through a characteristic velocity: 



M±v c . 



(3) 



If we consider the stell ar mass function fixed, we 
might take Af* ~ O.5M ; iMatzner & McKed (|l999al ) 
adopt v c ~ 40kms~ 1 in protostellar outflows, imply- 
ing 1 ~ io 39 - 6 



gems . 
Given M*, we may define 



(4) 



where SFRff is the star formation rate nor malized 
to the free fall rate at the mean density. iMcKed 
( 19891 ) estimates this parameter under the assump- 
tion that ambipolar diffusion determines the rate 
of collapse, and finds SFRff ~ 0.08 in regions 
well shielded from external far-ultraviolet radiation. 
The turbulence-regu l ated star formation theory of 
Krumholz k McKeel (|2005l ) predicts SFR ff ~ 0.034 
in nearby star-froming clumps. 

Although we used M* in defining v c and SFRff, it 
is important to realize that outflow-driven turbulence 
is invariant under the transformation M* — * aM*, 
SFRff — > oSFRfj, and v c — > a~ 1 v c , which leaves S 
and I unchanged. Therefore introduces no new 
dimensional scale so far as turbulence is concerned. A 
corollary is that, althoug h v r enters explicitly into the 
star formation efficiency ( Matzner Sz McKeel 200d ) . it 
will always appear in the combination SFRfff c oc SI 
in dimensionless ratios given here. 

Self-gravity and finite radius, on the other hand, in- 
troduce the dimensional quantities G (Newton's con- 
stant) and R, which can be used to form two dimen- 
sionless ratios with S, 1, and po. Useful forms include 
the radius ratio 

n* = t, (5) 



the acceleration ratio 



f 2 
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Po v 2 K 
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(6) 



which measures the acceleration scale SI/po in units 
of the surface gravity v 2 K /R = GM/R 2 (for Kepler 

l/2 n 7/2 



n: 



R 



speed Vr), and the quantity 2 esc /I 
which we define in § 14.21 

The scalings in equations ([I]) and pertain to an 
infinite (TIr 3> 1) and non-self-gravitating (II acc 3> 1) 
medium rather than a realistic clump or cloud. Con- 
sider the opposite: if Ur <C 1, then the region will be 
cleared by a typical outflow; and if Il acc <C 1, then 



outflows cannot provide the dynamical pressure re- 
quired for cloud support. One might therefore expect 
both parameters to be ^ 1 in a real cluster-forming 
environment. More specific predictions are provided 
in § 14.21 where we detail a model for turbulence in a 
cloud of finite radius and mass, and in §[7| where we 
consider the cloud's dynamical expansion or contrac- 
tion in response to outflow-generated turbulence. 

2.2. Finite sound speed and magnetization 

Continuing in the spirit of dimensional analysis, let 
us allow the gas to have a finite sound speed c s ; the 
relevant dimensionless parameter is 



IF 



Cot 



(7) 



which is an estimate for the inverse of the turbu- 
lent Mach number on scale I. Turbulence is super- 
sonic and strongly compressive so long as II C <C 1; 
for the parameters adopted above, this is true if 
SFRff > 10" 3 - 2 n^ 1 4 /6 (c s /0.2kms- 1 ) 7 / 3 , where the 
hydrogen density is 10 4 n#4 cm -3 . So long as II C <C 1 
we expect outflow-driven turbulence to be insensitive 
to the precise value of c s ; this is supported by the 
more detailed calculations given below in § 14.11 

The influence of sound speed on the star formation 
process itself should not be understated, however. By 
setting the thermal Jeans scale, finite sound speed 
deter mines both the normaliza tion of stellar masses 
(e.g., iPadoan fc Nor dhmdl I2002T ) and the star forma- 
tion rate (at least in the model of lKrumholz &: McKee 
2005), thereby affecting <S and X. 

The dynamics of the g as will certainly be affected 
by magnetic fields (e.g., McKee et al. 1993). These 
introduce at least two additional dimensionless pa- 
rameters: the degree of magnetization, as measured 
by the global mass-to-flux ratio or the local plasma-/? 
parameter; and the ambipolar diffusion rate, as mea- 
sured by the ratio of the neutral-ion collision time to 
the free-fall time, say. We expect that, like II C , these 
parameters influence the coefficients that appear be- 
low - the coupling factor A defined in §|H and the 
critical value of a va in §0- but introduce no other 
qualitative change to the theory. This expectation 
can only be tested by comparing the theory to obser- 
vations and simulations. 

3. OUTFLOW MOMENTUM PER STELLAR MASS, v c 

The outflow momentum per u nit stellar mass, v r , 
was es timated to be 40kms _1 by lMatzner &; McKee 
(200(3) on the basis that the typical protostellar wind 
velocity is ~ 200 km s -1 and about one-sixth of the 
accreted mass is ejected as wind - a value which is 
intermediate among several theoretical predictions. 
In order to check our estimate of v c , we combine 
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the Palla &: Stahler ( 19921 ) evolutionary models for 
ac creting protostars w ith an observational inference 
bv lRicher et al. (120001 1 of the relation between stellar 
and outflow properties. For those stars whose lumi- 
nosity is dominated by accretion, IRicher et al. show 
that the outflow momentum is roughly 0.3V}-*, where 
Vk± is the Kepler velocity at the surface of the prote- 
st ar. We obtai n Vk+ from iPalla fe Stahlerl 's models, 
and apply the IRicher et al.1 rule to all the material 
accreting onto a star of final mass M* in order to cal- 
culate its net outflow momentum 2~(M*). Dividing 
the total momentum by t he total mass f or a stellar 
population drawn from the Kroupal ( 200ll ) IMF with 



an upper cutoff of M* u , we find 



~ (28, 27, 22) log 



10 



0.11 M Q 



(8) 



assuming stars all accrete at (10 _5 ,3 x 
10 -5 , 10 -4 ) Mq yr _1 , respectively. The result is lower 
for a higher accretion rate, because these stars have 
less time to contract during accretion. 

For an upper cutoff M* u = 120 Mq this exercise 
gives about twice our f iducial estimate of 40 kms~ 1 , 
as previously noted by iTan &, McKeel ([2002) . How- 
ever, the excess momentum comes from massive stars 
that ( 1) require an extrapolation of the IRicher et al. 
(120001 ) scaling, (2) form only sporadically, relative 



to low-mass objects, and (3) may not form at all, 
if the cluster in question is too small to sample the 
IMF. For simplicity we hold to the initial estimate 
v c = 40kms~ 1 for the purpose of making estimates, 
although a more complicated model could be accom- 
modated using the theory we present in § 14.11 

4. VELOCITY SPECTRUM AND LINE- WIDTH-SIZE 
RELATION 

The winds that drive outflows possess a high energy 
per unit mass, which reflects the depth of the poten- 
tial well from which they were launched. Most of this 
is lost to radiation at the wind deceleration shock. 
What remains is the kinetic energy of the expanding 
shell, but this declines in proportion to the shell ve- 
locity as additional mass is swept up. At a radius of 
order I the shell loses coherence and merges with su- 
personic turbulent motions - which themselves dissi- 



pate e n ergy, as described, for i nstan ce. bvlStone et al.l 
(|l998T ). lMac Low et al.1 (|l99Sf ). and lMac Lowl (|l999T ). 
In contrast to the Kolmogorov cascade of incompress- 
ible turbulence, this supersonic cascade is not de- 
scribed by a constant energy flux to small scales. 

On the other hand, strongly radiative outflow shells 
conserve momentum in each direction, and therefore 
conserve a net scalar momentum X as they expand - 
where by "scalar momentum" we mean 



if p and u are the local density and fluid velocity, 
u cm is the center-of-mass velocity, and the integral 
extends over the volume affected by the outflow. A 
population of expanding outflows delivers momentum 
at a rate SI per unit volume. Turbulence is there- 
fore driven, on scale t, by a characteristic acceleration 
S1/pq. 

Does a radiative, supersonic cascade transport mo- 
mentum to small scales the way an incompressive 
cascade transports energy? This seems plausible, 
based on the dynamics of expanding outflow shells. 
Moreover, this hypothesis is consistent with the well- 
known scaling laws of supersonic turbulence. The 
one-dimensional (spherically averaged) spectrum of 
the velocity field u(x) is defined as 



l r°° 

E ( k ) = 7T" / ( u ( x ) • u ( x + nr)) e 
27T Jo 



ikr 



dr (9) 



where angle brackets indicate an average over time, 
space, and the direction of the unit vector n. Sup- 
pose the turbulence transports an acceleration SI/ po 
from £ to small scales. As S has units cm -3 s -1 and 
as ST / po has units cms -2 , 



CT 

E{k)=A—k- 2 
Po 



(10) 



for some coefficient A which depends on other di- 
mensionless parameters. Equation (|10p matches the 
spectral slope for supersonic tu rbulence, as deter - 
mined by numerical experiment ( Porter et al. 19921 ). 
It coi ncides with the spectral slope in Burgers' turbu- 
lence (lKidal|l979T ). and indeed with that of individual 
shocks. Moreover it resembles lLarson 's (119811 ) and 
matches Sol omon et al. I's (119871 ) estimates of the line- 
width-size relation in molecular clouds. Therefore, we 
adopt the central hypothesis that outflow- driven tur- 
bulence transports the flux of scalar momentum per 
unit mass (i.e, an acceleration) from the driving scale 
to small scales, where it is destroyed by the collision 
of oppositely-directed motions. The remainder of this 
paper amounts to an exploration of this hypothesis. 

Consider the one-dimensional velocity dispersion, 
or line-width, a(r). To be specific we define 3<r(r) 2 /2 
to be the mean kinetic energy per unit mass of spheres 
of radius r, relative to their centers of mass. Since for 
now we consider the gas infinitely cold (an assump- 
tion we drop in §[3J), a — > as r — > 0. On small scales 
(r <C £), o~(r) is determined by dimensional analysis: 



p\u 



u, 



|d 3 x 



o-(r) 2 = A 2 — r (r < 
Po 



for some coupling coefficient A. 



(11) 
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A lower limit on A can be estimated from the fact 
that o~(r) 2 includes the kinetic energy of decelerat- 
ing shells. The mean kinetic energy density in shells 
traveling faster than v is 



„ /■*(") Iv' , , 

S / dt' 

Jo 2 



(12) 



where t'{v') is the age at which the flow decelerates to 
v'. Since the shell expands according to dr' = v'dt', 
the integral can be written 



r(v) 



dr' 



r(v) 



SI, 



(13) 



which is SIi/2 at the merging scale, r(v) = r{a) = i. 
The total energy density on this scale is 3a 2 (£)/2; 
comparing to equation ([TT]). we find 



A 2 > 1/3. 



(14) 



This is only approximate, since our estimate of the 
kinetic energy did not account for the complicated 
dynamics of merging, and since equation (jlip is ex- 
pected to break down at the merging scale. However, 
it guides us to expect A 2 to be of order unity, and to 
attribute a value in excess of 1 /3 to the persist e nce o f 
turbulent energy. In fact, the lKrumholz et al. I (|200d ) 
estimate of turbulent driving and decay corresponds 
to a change in A 2 of 0.64 for this persistent turbu- 
lence, so our best estimate of A 2 is 1/3 + 0.64 = 0.97. 

On scales large enough that no individual outflow 
can dominate the motion, the velocity spectrum is 
likely to be affected by a cascade from even larger 
scales. Indeed, if outflows are sufficiently weak or the 
driving is sufficiently strong, it is possible for this cas- 
cade to dominate over outflow-driven motions even on 
scales smaller than i. In the equations below we shall 
include a term to describe this external cascade. 

4.1. Variations of outflow intensity 

Up to this point we have considered outflows of uni- 
form intensity. However fluctuations should be com- 
monplace - due to variations in the mass and accre- 
tion rate of the driving star, for instance, or due to 
the angle from which a collimated outflow is viewed. 

To treat this possibility, consider a single outflow 
whose momentum scale I\ differs from the character- 
istic value I. It dominates the motion of gas it over- 
takes within a "merging" radius r m i, defined such 
that 

M(r ml )a(r ml )=I 1 . (15) 

This equation relies on the basic assumption that the 
one-dimensional velocity dispersion a(r) can be used 
to gauge the merging radius, without reference to the 



intermittency of the motions on scale r. Moreover it 
neglects magnetic forces, which can cause merging at 
smaller radii. Here M(r m \) refers to the mass en- 
countered within a radius r m \ of the driving source; 
we write 

Mir) = Y<PmPor 3 (16) 

where 4> m > 1 if stars tend to form within extended 
overdensities. In general <p m is a function of the other 
dimensionless parameters, including r/£, but we shall 
treat it clS 3j number. 

We previously considered the outflow strength I 
and rate S as unique values. We generalize now to a 
range of values for I, and consider <S(< I) as the cu- 
mulative rate of outflows weaker than I. To proceed 
we suppose that a(r) 2 - which is proportional to the 
specific kinetic energy of motions on scales smaller 
than r - is a linear superposition of the contributions 
from outflows of a single intensity. 

Therefore, let us reexamine the single-intensity 
limit considered in §[2j This is given by S(< I) — > 
S\H(I — Zi), where H(x) is the Heaviside function, 
and where subscript '1' denotes a single value of Z; 
then dS/dl = S\5(I — Ii). The contribution o~ 2 (r) 
is oc r on scales smaller than r m i, and is constant on 
larger scales; therefore we adopt the idealization 



a 2 (r) = A 
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r < r m i 
r > r m i 



(17) 



The contribution of outflows to a\ on scales larger 
than r m \ should not be interpreted as an inverse cas- 
cade: it simply reflects the fact that o 2 (r) incorpo- 
rates all motions within spheres of size r. Differenti- 
ating twice, 



d 2 a\ 
dr 2 



A 2 S 1 I 1 
Po 



S(r - r ml ) 



A 2 Zi dS/dl 
Po dr ml /dl' 



This equation holds for each value of Z in our linear 
superposition model, so in general 



2^2 



d z a 
dr 2 



A 2 I dS 
Po dr ' 



(19) 



In this equation, as in the rest of this section, we 
take Z to mean the value for which r m \ = r according 
to equation (|15p . so that Z, 5, and r are monotonic 
functions of one another; dS/dr could be replaced 
with [dS / dl){dl / dr). Keep in mind that the rela- 
tions between these quantities depend on <r(r), which 
must be determined self-consistently. 

Integrating, 



da 2 A 2 f Sto " , IO/ 

— — = a ex t / IdS 

dr po Js(r) 



(20) 
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where a ext is the acceleration due to a cascade, if any, 
from larger scales, and «St t * s ^ ne t°t a l outflow rate 
per unit volume. Therefore 



, da 2 
lim — — 

r— >o dr 



flext + 



A 2 
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IdS, 



(21) 



which generalizes equation (jlip . Taking a 2 to include 
thermal motions, we also have 



lim a 2 

r^O 



(22) 



Given a total momentum injection rate JXdS, 
what form for <S(< X) maximizes a(r) on large 
scales? It is obvious from equation (|2j) that since 
i/t oc (IS) 3 / 7 ! 1 / 7 , turbulence is enhanced if momen- 
tum is emitted in fewer, stronger bursts. (Up to a 
point: in a finite region, bursts that are too strong 
will escape, and those that are too rare are irrele- 
vant.) Suppose we fix the number of events and thus 
the mean outflow strength 



X 



J IdS 



'tot 



as well as / XdS. Then, by the same argument, tur- 
bulence is maximized if all but one outflow is negli- 
gible in strength, and the last carries all the momen- 
tum. 

More generally, we expect that turbulence is en- 
hanced by broadening the range of outflow strengths. 
Consider the integration by parts of equation (|20p 
over radius, from to r: 



a 2 (r) 



+- 



A 2 

+ ra ext H r 

Po 

a 2 r s ( r ) 
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S(r) 
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r'l'dS'. 



(24) 



This solves differential equation (|19p , but only implic- 
itly, since the function r(X) given by equation (|15p . 
and therefore <S(r), depend on a(r). Equation (f24"|) is 
useful, however, because the first integral disappears 
at the maximum outflow extent (r = maxr m i). The 
contribution of outflows to a 2 (r) on even larger scales 
is therefore given by the second integral 



rldS. 



(25) 



A 2 
P Jo 

This allows us to define a turbulence enhancement 
factor as the ratio of the outflow-driven turbulent 
pressure (eq. [25] 



to what it would be for a mono- 
lithic model in which X = X: 



/ 5tot rXdS 



£ 



47T0„ 
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(26) 



This factor illustrates the increase in a 2 expected due 
to the broadening of the outflow distribution. It de- 
pends on c s and a ext because of the implicit nature 
of equation (|24p . but to be definite we calculate it 
assuming c s = a ex t = 0. 

Our expectation that broader distributions of out- 
flow strengths lead to stronger turbulent motions is 
verified in Table [H where £ is presented for the out- 
flow distributions discussed below. 

4.1.1. Mass function 



We now consider the two primary sources of vari- 
ation in outflow intensity: differences in the stellar 
mass, in this subsection, and differences with angle 
in collimated winds in the next subsection. 

Let F(> M+/M+) be the fraction of stars that 
exceed the mean stellar mass M± at birth; then 
-F(O) = 1 and J x\dF/dx\dx = 1 if the limits of in- 
tegration include all stellar masses. Assuming X is 
strictly proportional to M+ as in equation ([3|), and 
parametrizing the star formation rate using equation 
(23) W e have 



3W 



M*v c 



(27) 

The corresponding <r(r) 2 achieves a value enhanced 
by £ = 1.85 relative to a single-mass model. 

4.1.2. Outflow collimation 

A second and even more significant source of vari- 
ation in outfl ow intensity is the collimat ion of proto- 
stellar winds. Matzner Sz McKeel ( 1999al ) argued that 
magnetic stresses cause both disk winds and X-winds 
to approach the same asymptotic structure at large 
distances from the source: the apparent strength at 
an angle p to the outflow axis is 

z( m )=xp(H), p(H) 



ln(2/0 o )(l + o 2 -^ 2 ) 
(28) 

where the normalization factor ln(2/#o) ensures 
/ P(|/i|)<i|/x| = 1 so long as #o <^ 1- The soften- 
ing angle #o encompasses precession, internal shocks, 
fluid and magnetic instabilities, phenomena at the 
jet-ambient interact ion, and anything e lse tha t de- 
grades collimation. Matzner &; McKed (|l999aT l esti- 



mate #o ~ 10 for several well-observed outflows. 
This implies a variation of ~ 10 4 in outflow strength 
with angle, which has a strong effect on the ampli- 
tude and character of outflow-driven turbulence. (It 
is possible, however, that a larger value of #o applies 
to the later stage of outflow evolution and to the driv- 
ing of turbulence.) 
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Given an outflow of net momentum Z, the apparent 
strength equals I at angle 



P-\X/X) 



1 + 



ln(2/0 o )X/X. 



1/2 



(29) 

where P 1 is the functional inverse of P, so long as 
this gives < \fj,\ < 1. For a random orientation 
is evenly distributed within this range; therefore 
P~ l (I/I) is the cumulative probability that the ap- 
parent strength is weaker than 2", and dP -1 / dil /I) 
is the distribution of I/I. (Formally, dP" 1 /d(I/I) = 
d\fi\/dP.) 

We adopt the basic assumption that turbulence 
driven by collimated outflows is very similar to tur- 
bulence driven by spherical outflows, if one general- 
izes the distribution S(I) to the apparent distribution 
S(I), which is the convolution of iS and P _1 over the 
logarithm of their arguments: 



dS 



dS 



dP- 



dlogl dlogl d\og(l/I 



(30) 



The influence of collimation on a{r) is illustrated 
in figure [H which displays numerical integrations of 
equation (|20|) . Two points can be drawn from these 
results: (1) A finite thermal sound speed is of no prac- 
tical consequence so long as i > c s t, since solutions 
converge when a > c s ; (2) a(r) continues to rise until 
r ~ max(r m i). (Recall that max(r m i) refers to the 
value obtained by setting I — > maxl.) However the 
rise in <r(r) slows to zero as this limit is approached. 

We expect the coupling coefficient A to depend 
somewhat on the degree of collimation, thanks to 
complicat ed phenomena like c o llision s between young 
outflows (jCunningham et all 120061 ). although the 
sense and magnitude of this dependence are un- 
known. 

4.2. Finite-size regions 

What if the region of interest (a dense "clump", 
say) is too small to catch all the outflow momentum? 
The eruption of outflows from protocluster regions is 
frequently observed, and this p r ocess deserves special 
attention. IMatzner fc McKed {2000) calculated the 
mass ejection rate in this case; we wish to consider 
the outflows' dynamical effects. 

It is not sufficient to simply evaluate <r(r) at the 
clump size R using the results of the previous section. 
Equation (i20j) assumes a homogeneous background 
on scales larger than r, as a(r) reflects the downward 
cascade driven by outflows merging on larger scales, 
as well as the upward cascade composed by those 



outflows as they expand. However, outflows take mo- 
mentum as well as mass with them when they escape. 
Collimation allows this to happen in some directions 
without the entire clump being unbound. Specifi- 
cally, escape occurs in direction \i when X(|/x|) > 2"eso 
where 

Xesc = c g Mv esc (31) 

if M = M{R) and v esc = (2GM/R) 1 / 2 are the mass 
and escape velocity of the region. The factor c g ac- 
counts for gravitational deceleration of an expand- 
ing shell, which saps its momentum. It is close to 
unity, however: for densit y distributions p(r) oc r~ kp , 
Matzner fc McKed (l2000l ) found 

1/2 



9-3^ 
8 - Zk, 



(32) 



(their equation [A13]). This formula holds when out- 
flows are driven impulsively - a safe assumption, so 
long as the wind that drives an outflow has a du- 
ration similar to the free-fall time of its collapsing, 
overdense core. 

Given that outflow intensities exceeding I esc get 
away, the theory of § 14.11 is easy to modify: simply 
replace <St t with S(I esc ), throwing away the remain- 
ing momentum. It is important to realize that some 
outflows emerge from the clump surface and rain 
back on it later; the replacement just suggested treats 
them no differently from those that merge within the 
clump. Although approximate, this approach cap- 
tures the essential division between capture and es- 
cape. 

The effect of outflow eruptions on the velocity scale 
o~(R) is determined by how much outflow momentum 
is eliminated in this procedure. It depends, there- 
fore, on the dimensionless ratio I esc /I as well as on 
the shape of S(I) (see also § 12.11) . This dependence 
is illustrated in figure [2j where a slow dependence on 
I esc /I is apparent. This is natural, because the ap- 
parent intensity I can exceed T by a factor of 10 6,7 
in the model plotted - of which 10 2 ' 7 comes from the 
range of stellar masses, and 10 4 arises from collima- 
tion. F o r a d efault model in which we adopt the 
Kroupal (j200lh IMF, take v c to be independent of 



M*, consider winds to be collimated with 9q = 10 



and take a. 



cxt 



0, we find 



a 2 t 2 



c 2 t 2 



+ A 



+ K, 



-277 



-l/r, 



(33) 



where the fit parameters K\ , , and r? are functions 
of Iesc/I: 

-1/10 



2.33" 10 + ( 0.32 + 0.40 log 



10 



I 



-10' 



(34) 



8 



2.5 



1.5 



< 
b 



0.5 



10 ' 



Kroupa IMF,e Q =10~ 



Kroupa IMF, spherical 



10 



10" 

r/l 



10 



F ig. 1. — - Line- w idth er(r) versus size, normalized to the characteristic values in § for outflows driven by stars drawn from 
the iKroupal (|2001f ) IMF. If outflows are spherical, the theory of 5 14.1. II gives the lower curve; if they are collimated, the theory 
of S I4.1.2I gives the upper cur ve. At scales much smaller t han the driving scale £ , a oc r 1 ^ 2 - in agreement with the Larson's-law 
scalings for molecular clouds (|Larsonlll98lD as reported bv lSolomon et aO ll 19871 ). 



Table 1 

Turbulence enhancement factors. 



Mass function 


£ (isotropic) 


£ (9 = lO" 2 ) 


Single-mass 


1.0 


3.04 


KrouDa f 20011 IMF 


1.85 


5.38 
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K 2 



and 



1.22-°+ 0.55 + 0.26 log 10 



T 



-6 



-1/6 



ll + 151og 10 (X esc /X) 



(35) 
(36) 



l + 101 Ogl0 (XescA) ' 

This fit reproduces the numerical evaluation of a(r) 
to within 5% for I esc > 61. The quantities used in 
these formulae are, in convenient form, 



fc = 770 Ml'^^^-, 



1/7 



and 



(37) 

pc, 
(38) 



^ = 1.24M- 1/28 S1 28 (- 



4 / 7 / 0.034 \ 3 / 7 km 



t 3 " cgs Uokms- 1 ^ VSFR ff y s 

(39) 

where M(R) = 10 3 M 3 M Q and M(R)/(irR 2 ) = 
E cgs gem -2 , and we used k p = 1.5 to evaluate c g . 

5. COMPARISON TO SIMULATION 

The simulation bv lLi Nakamural {2006) offers the 
closest comparison to the theory presented above. 
These authors simulate a self-gravitating, weakly 
magnetized, initially turbulent cloud, identifying re- 
gions of localized collapse and replacing them with 
spherical outflows. Using their fiducial conver- 
sion from code units to physical units, their out- 
flow momentum scale is nearly constant at X = 
10 39 ' 7 gcms _1 . The 960 M & of cloud gas is initially 
arranged in a mildly centrally peaked density pro- 
file of outer radius 0.75 pc, and becomes more cen- 
trally condensed over the course of the run. After 0.6 
Myr star formation begins in earnest, and between 0.9 
and 1.2 Myr the stellar fracti on increases from 6.5% 
to 13.4% of the cloud mass. iLi k Nakamura report 
a final, three-dimensional rms velocity of 1.5 kms -1 
corresponding to a = 0.85 km s -1 . 

If one assumes that the gas has contracted by a 
fraction fn of its initial radius during the period 
0.9-1.2 Myr, then p ~ 10~ 19 ' 4 / R 3 g cm -3 and S ~ 
10 



-66.6 



f R cm 3 s in that period. Equation (|17p 
then predicts a = 1.7/RAkms -1 , which is consis- 
tent with the lLi k Nakamural results if A = 0.6f R 3 ^ 7 . 
Contraction by fn = 1/2 appears reasonable based 
on their figure 2, in which case 

A ~ 0.8. 

This seems entirely consistent with our expectations 
from §[2] that A ~ 1. Note, however, that one should 



expect changes in numerical resolution, cloud magne- 
tization, and outflow collimation to be accompanied 
b y variations in A . 

ILi k Nakamura also mention that collimated out- 
flows appear to drive more vigorous turbulence, a s 
we would expect from § 14.1.21 Nakamura k Lil (2006) 
have modified their previous simulations by collimat- 
ing outflows within 30° half opening angle, and report 
that the enhancement is not discernible from the sim- 
ulations. In our theory, this modification amounts to 
increasing X by a factor of 7.5 while holding ST fixed; 
we would expect a to increase by 7.5 1 / 7 = 1.3 as 
in eq. [2]), so long as A does not change. Whether 
this mild enhancement of turbulent velocity is seen 
in simulation will require further investigation. In 
an equilibrium cloud, as our referee points out, one 
must tease this effect from the redistribution of mass 
caused by more vigorous turbulence. 

A deeper i nvestigation of t h e sam e class of model is 
reported by iNakamura k Lil ( 20071 ). who report sev- 
eral results in line with our predictions. In this work, 
outflows posses a 30° "jet" component and a uniform 
"wind" component. The authors report that the jet 
component is more effective at driving turbulence, 
based on the slowdown of star formation observed 
if jets are strengthened (all else being equal). Their 
turbulent power spectrum obeys E(k) cx k~ 2 as antic- 



ipated in equation (TTOj) . and flattens at large scales - 
possibly reflecting the flattening of a{r) (figure [1]). 
Finally, they obtain an equilibrium density profile 
p cx r~ L5 , in line with the argument we present below 
i n m 



Mac Low! (|2000h has also simulated turbulence 



driven by collimated sources, but a direct compar- 
ison is not possible because the driving field in his 
simulation was steady rather than impulsive. 

6. COMPARISON TO OBSERVATION: NGC 1333 

Located within the Perseus molecular cloud, the 
NGC 1333 reflection nebula is the site of a vigor- 
ous burst of clustered star formation within a dense 
molecular clump. The current population of at least 
143 young stars is estimated to be 1-2 Myr old 
(|Lada et alJ [l996K and about ten molecular outflows 
and Herbig-Haro syste ms emerge from the most re- 
cently formed objects. lLada k Lada ( 20031 ) estimate 
a stellar mass of ~ 79 M©, based on the limiting K 
magnitude of 14.5. Distance estimates range from 
D =210 pc to 350 pc, as discussed by IWarin et al. 



(1996). 



Eldge et all (|2003l ) map the NGC 1333 core (and 



many others) in 13 CO(1-0), C 18 O(l-0), and C 18 0(2- 
1), probing sequentially smaller and denser regions. 
They derive masses from a large velocity gradient 
analysis; combining these with their reported veloc- 
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F ig. 2. — Effec t of finite clump radius on turbulence driven by collimated outflows (do = 10 2 ) driven by stars drawn from 
the iKroupal 1)200 IT ) IMF. The thick curve represents an infinite uniform medium; the other curves are labeled by I csc /X, which 
determines how much momentum is lost from a finite clump. Circles mark the edge of the region for fiducial parameters uhi — 1, 
v c = 40 kms -1 , SFR ff = 0.034, and A = 1. 



ity dispersions, and assuming that gas dominates 
the mass budget, we d erive the virial parameter 
(|Bertoldi fc McKedll992l ) 



a. 



5Ra(R) 2 
GM(R) 



[1.27, 1.32, 1.19) (40) 



for ld CO(l-0), C ls O(l-0), and C ls O(2-l), respec- 
tively. Since a v ; r is thought to vary, in strongly 
self-gra yitating regions, betw een 1.11 (for molecular 
clouds; ISolomon et al] Il987h and L34 (for massive, 
magnetized molecular cores; McKee fc Tanll200al ). 
adopt D = 220 pc (close to D = 212 pc adopted by 
Quillen et all 120051 ). The embedded star cluster then 
have an effective radius R = .39 pc, scaling the value 
given bv lLada fc Lad3 (|2003l ) to the closer distance; 
the enclosed mass is then 446 M , and the escape ve- 
locity is 3.1 kms -1 . The velocity dispersion at the 
same radius is a = 1.1 km s -1 . 

Can this velocity dispersion be maintained by the 
protocluster outflows? To eval uate the formu lae in 
§0 and g21 we adopt the IKroupal (|200ll ) IMF 
and assume the current star formation rate equals 
79 Mq per Myr. Other model parameters are <Stot = 
10- 65 - 8 cm- 3 s _1 , po = 10 
10" 2 , and X PSC = 10 415 



-18.9 



gems 



gciir^ T = 20 K, 9 = 
1 (where we have used 



k p = 1.35 to estimate c g = 1.1 in eq. [32]). For ref- 
erence, these values give SFRff ~ 0.034 in the region 
of interest - virtually id entical to the value SFRff = 
0.038 we derive from the iKrumholz fc McKed (feOOSl ) 
theory. We take a ex t = to probe the influence of 
outflows alone. 

Using v c = 40kms _1 and applying the prescrip- 
tion of §4.21 to predict the velocity dispersion at 
R = 0.39pc, we find a(R) = 1.12(A/0.8) kms -1 . 



This is entirely consistent with the observed value 
of 1.1 kms -1 , for precisely the same value of A we es- 
timated in §0 This result is entirely consistent with 
the proposition that turbulence in the star-forming 
region of NGC 1333 has been regener ated by out- 
flows. The cluster outflow inferred by IWarin et al 



(1996) could then legitimately be viewed as the mass 
ejected by e r uptin g jets, and the inflow detected by 
I Walsh et al.1 ( 20061 ) might represent a return flow of 
gas ejected below the escape speed. 

Although the value v c = 40 kms -1 (which was 
justified in § [3] using the ob servational relations re- 
ported bv lEicher etaT1l2000l ) is in line with the mo- 
mentum of the NGC 1333 outflow HH 7-11 as es- 
timated bvl Snell fc Edwards! (Il98lh . iKnee fc Sandelll 
(|2000r ) and iQuillen et al.1 (|2005l V have inferred signif- 



icantly lower m omenta for the NGC 1333 outflows. 
Knee Sz Sandelll estimate ~ 1M km s - 1 per outflow 



for ten currently active outflows, and IQuillen et al. 



derive a similar momentum scale for the explosions 
that open cavities i n the ga s . Giy en that the mean 
stellar mass in the IKroupal (|200ll ) IMF is 0.21 M Q , 
this momentum scale corresponds to v c ~ 5 kms -1 . 
In the formalism of § 14.21 this lower value would imply 
a(R) = 0.47Akms -1 . 

It is possible that, for strongly collimated out- 
flows in a real, magnetized star-forming region like 
NGC 1333, A = 2.3 as required to make this 
lower estimate consistent with the observed value of 
1.1 kms -1 . However it is equally probable that the 
inference v c ~ 5 kms -1 is an underestimate: com- 
pare this to the estimates in §[3l and note that the 
minimum value estimated by iRicher et al. (2000, m 
their SIII.B) is ~30kms _1 . One solution to this dis- 



11 



crepancy could be that the outflows now observed in 
NGC 1333 are simply smaller than the avera ge - pos 



sible, since the median value of M* in the iKroupa 
(|200lh IMF is 2.6 times lower than the mean value. 
Alter natively these outflows could be more powerful 
than iKnee k Sandelll (|2000l ) report. The fact that 
they find a 2.5M0kms _1 as the net momentum of 
HH 7-11 in their 12 CO(3 - 2) an d 12 CO(2-l) survey, 
whereas Snell &: Edwards ( 1 98 ll ) find ~ 34kms _1 for 
the same outflow in observations involving the four 
transitions 12 ' 13 CO(J=2-1,1-0), raises this possibility. 
CO optical depth, uncorrected inclination, motions 
at velocities close to systemic, and loss of momen- 
tum from the cloud can all c ause outflow momentum 
to be underestimated. (See IWalawender et aL 2005 
for a discussion of some problems in the determi- 
nation of outflow momentum.) Moreover, the most 
noticeable outflows are the young, rapidly expand- 
ing ones. Since these are still being driven, they will 
not have acquired their final momenta. Finally, note 
that very little mass can be ejected from the region 
by outflows if v r is as low as 5kms _1 : the theory of 
iMatzner fc McKed (l2000h predicts a star formation 
efficiency e ~ 88% in this case, whereas e — 47% if 
v c = 40kms _1 . The latter value is more consistent 
with the fact that only ~ 18% of the mass is currently 
in stars. 

For all of these reasons we consider it most likely 
that turbulence in NGC 1333 is driven by outflows 
with v c ~ 40kms~ 1 (such that A ~ 0.8 in the theory 
of § I4.2p . although the alternatives - that A ~ 2.3, 
or that outflows are only a minor contributor to the 
turbulence - cannot be ruled out. The last option 
is especially unattractive, as it leaves unanswered the 
question of how turbulent energy is regenerated there. 

7. IMPLICATIONS FOR STAR CLUSTER FORMATION 

Our analysis of NGC 1333 supports the assertion 
that outflow-driven turbulence affects the dynamical 
state of gas clumps while star clusters form within 
them. How, then, does feedback influence the cre- 
ation of star clusters? 

A preliminary answer can be gleaned from figure 
O where we calculate a v ; r for cluster-forming re- 
gions of various total mass (M) and column den- 



sity (E = M/{%R)) using the formalism of § 14.21 
In this calculation the region was assumed to be 
neither expanding nor contracting, and devoid of 
any additional source of turbulence (a ex t = 0); out- 
flows were collim ated with 9n = 10~ 2 and v c = 
40kms~\ and the lKroupal (l200ll ) IMF was adopted. 
In the left panel, the star formation rate parameter 
SFRff was held fixed at the value 0.034 we derived 
for NGC 1333; on the right, SFRff and a v f r were 
derived self-consistently according to the model of 



Krumholz fc McKed (|2005h : however the gas temper- 
ature was held fixed at 20 K for simplicity. 

We expect that there exists a critical value of a v i r , 
between 1 and 2, such that a cluster- forming clump 
will be in virial equilibriu m given a mode r ate e xter- 
nal pressure; for instance, iMcKee fc Tail (l2003h de- 
rive a v ; r = 1.34 for their equilibrium model of mag- 
netized, turbulent cores. Smaller values of a v ; r imply 
contraction. Higher values imply either expansion, or 
require a strong confining pressure. 

An immediate conclusion from figure [3] is the exis- 
tence of an outflow-driven equilibrium state (a v ir — 
1.5) for column densities and masses relevant to star 
cluster formation. It is however an unstable equi- 
librium: since 3a v ; r /9S < at fixed mass, turbu- 
lent support is weakened if an equilibrium clump is 
compressed. The instability is reduced - though not 
removed - by the self-regulation of the s tar forma- 
tion rate in the iKrumholz & McKed (120051 ) model. It 
would be reduced even further if we were to account 
for heating of the gas by star for mation, as the sound 
speed enters (weakly) into the IKrumholz &: McKee 
theory. 

Unstable equilibria were previously found by 
iMcKed (1989) in models for gian t molecular clouds, 
and_WnMatzner| (jl999l ) and IMatzner & McKed 
(|l999bl i in simpler models for star cluster formation. 
In fact, instability is generic to any model in which 
SFRff and v c are slowly varying functions of the cloud 
parameters. For a virialized cloud this can be seen by 
comparing the cloud dissipation rate ~ a(R) 5 /G with 
its energy input rate - which scales as the star forma- 
tion rate ~ SFRff a(R) 3 /G, times the energy injection 
per unit star mass ~ v c a(R). So long as SFRff v c does 
not increase as fast as a 2 when the cloud contracts, 
dissipation overcomes energy injection. Sharp thresh- 
olds, such as the photoionization column in IMcKed 
(1989), are therefore required for stable equilibria. 

The scenario for star cluster formation implied by 
figure [3] is qualitative l y cons istent with the one ad- 
vance d by iTan et aL I (|2006l ) and IKrumholz fc Tanl 
(2006]), in that turbulence supplied by outflows is 
important enough to slow collapse over several, per- 
haps many, crossing times. Because the equilibrium 
state is unstable , is al so consistent with proposal by 
Li fc Nakamural EH) that this support may fail in 



some circumstances, and that this might lead to a 
state of collapse that triggers massive star formation. 
However, massive-star outflows are almost certainly 
more powerful than our si mple model has assum ed - 
as discussed in § [3] and by ITan & McKed" (j2002f ). If 
there is a column density threshold for massive star 
formation, then it is possible that the feedback from 
these stars creates a stable equilibrium - in which 
E remains near the threshold. Alternatively, massive 
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stars might disrupt the clump outright; we leave these 
questions for future investigation. 

The form of the line-width-size relation cr(r) affects 
the internal structure and dynamics of a clump, in 
addition to its global expansion or contraction. Sup- 
pose we define (3{r) = din a/ din r as the logarith- 
mic slope of this rela tion. In a singular polytropic 
model for the clump ()McKee k Tanll2003l ). the den- 
sity index k p (in p oc r~ kp ) is related to (3 through 
k p = 2 — 2(3 (and both (3 and k p are assumed con- 
stant in radius). In our model, (3{R) can take any 
value between and 1/2, depending on the ratio be- 
tween R and the turnover in a{r) - as seen in figure 
[2j This figure also shows that star-forming regions 
with uh ~ 10 4 cm -3 have radii near this turnover, 
so that (3{R) ~ 1/4 - generally consistent with the 
value j3 = 0.21 ± 0.03 obtained fo r massive cores and 
clumps bv lCaselli k Mversl ( 19951 ). More specifically, 
the fit given in equation (|33p implies 



1 



2H 1 + 



cjt 2 
7JM 



H 



(41) 



H 



,'K 2 



(42) 



where 

To the extent that the IMcKee k TanI (|2003T ) model 
can be generalized to our model for turbulence, (3 ~ 
1/4 implies k p ~ 1.5 on the scale of the region. 
This e stimate coincides with the value McKee k TanI 
(|2003l 1 adopted after a survey of the observational lit- 
erature. If our proposal for the origin of k p is correct, 
then one would expect that k p is steeper (closer to 
2) in denser and more massive regions (so long as 
turbulence is dominated by outflows), as both R/l 
is higher there , and as I/2 esc is lower. Note that 



Huff k Stahieil (j2006h derive k p ~ 2 for the birth re- 



gion of the Orion Nebula Cluster. 

The line-width-size index (3 also enters into es- 
timates for the ty p ical r adii of protostellar disks. 
Kratter k Matznerl ( 20061 ) calculate the specific an- 
gular momenta of cores prior to their collapse, and 
find that a larger value of (3 leads to a greater an- 
gular momentum. Note that, since (3 is a decreasing 
function of r, its value i s greater on the core scale 
than on the clump scale. Tan et ail ( 20061 ) have pre- 
viously proposed that (3 ~ 1/2 on small scales, falling 
to ~ 1/4 on scales large enough to sample the tidal 
field. We propose that this shift is produced by the 
dynamics of driving - and shapes the gravitational 
potential, rather than reflecting it. 

8. CAVEATS 

The theory presented here relies on several assump- 
tions and approximations. The central hypothesis, as 



stated in $H is that a strongly dissipative, supersonic 
turbulence can be described as a cascade of (scalar) 
momentum to small scales in much the same way that 
Kolmogorov turbulence can be described as an energy 
cascade. This is not too bold, as it agrees with the 
spectral slopes observed in simulations of these cas- 
cades. In § 14.11 we assumed that these cascades will 
simply add, in the case that turbulence is driven by 
more than one type of source. In § 14.1.21 we assumed, 
further, that this superposition can be applied to col- 
limated outflows - if the variation of strength with 
angle is treated in the same manner as a population 
of isotropic sources of different strengths. Moreover, 
our treatment in § 14. 21 assumes that the turbulent ve- 
locity in a region of finite radius can be estimated 
by (1) discarding momentum ejected in escaping out- 
flows and (2) evaluating the resulting line-width-size 
relation at the scale of the region. Finally, we have 
assumed that the differences between outflow-driven 
turbulence in magnetized and unmagnetized gas can 
be encompassed by modifying the coupling coefficient 
A, and in §[71 by a modifying the critical virial param- 
eter corresponding to collapse. All of these assump- 
tions are sure to be approximate at some level, and 
must be tested against numerical simulations and ob- 
servations more thoroughly than we have done in §[5] 
and §[6j 

9. SUMMARY AND CONCLUSIONS 

Based on our observation that the spectral slope of 
supersonic turbulence implies a momentum cascade 
rather than an energy cascade, we have constructed 
a simple model for turbulence stirred by momentum 
injection from stellar outflows. A key feature of our 
model is that turbulence is enhanced if momentum 
injection is spread amongst outflows of a wide range 
of strengths, or, by extension, if outflows are strongly 
collimated. These effects allow for some of the mo- 
mentum to be deposited on relatively large scales 
where turbulent decay is slow. The diversity of out- 
flow strengths (reflecting the range of stellar masses), 
and their strong collimation by magnetic fields, both 
imply that real outflows are quite effective at driving 
turbulence. 

Our comparison to NGC 1333 supports the asser- 
tion that turbulence in the cluster-forming gas has 
been entirely regenerated by outflows. Although un- 
certainties regarding the outflow momentum scale 
and the value of the theoretical outflow-turbulence 
coupling coefficient A limit the strength of this con- 
clusion, these uncertainties will be eliminated by fu- 
ture numerical and observational studies. 

In our fiducial model, outflows will maintain tur- 
bulence within cluster-forming regions of a few thou- 
sand solar masses, whose column densities are ~ 



13 




'1 234561 23456 

log lo M(M ) Iog 10 M(M @ ) 



Fig. 3. — Virial parameter a v ii for outflow-driven turbulence within clouds of various total mass M and column density 
S = M/(nR 2 ) . Outflows are assumed to be collimated ((9o = 10 2 ) and driven, with v c = 40kms , by stars drawn from the 
iKroupal (|200ll ) IMF; A = 1. In fa), SFR ff is held consta nt at 0.034; in (b), SFR ft is derived self-consistently from civir using the 
star formation rate model of lKrumholz fc McKed ((2005) . In both cases the gas temperature is assumed to be 20K throughout, 
and there is no external driving (a ox t = 0). Conditions for expansion and contraction are estimated assuming that virial balance 
corresponds to a value of a v ir between 1 and 2. 



0.3 — lgcm -2 . Intriguingly, the energetic equilib- 
rium in these models is an unstable one. It is possi- 
ble that this leads to g lobal contraction or collapse, as 
iLi fc Nakaniui^l (j200fih have suggested. However this 
instability could equally be an artifact of our assump- 
tion that outflow momentum scales in proportion to 
the mass of the forming star. 

The structure of a turbulence-supported region re- 
flects its turbulent spectral slope. The radial density 
index takes values close to -2 if turbulence is driven 
on small scales, or close to -1 if it is driven on scales 
larger than the region of interest. Our model predicts 
intermediate slopes (~ —1.5) because outflow colli- 
mation allows both of these to occur simultaneously. 
We expect this slope, which is consistent with obser- 
vations of cluster-forming regions, to steepen in re- 



gions that catch outflow momentum more effectively. 

Although our model was developed to address pro- 
tostellar outflows, it could in principle be extended 
to other forms of driven, supersonic turbulence in 
radiative gas, such a s HH-reg i on-dr iven turbulence 
in molecular clouds (jMatzner 20021 ) or supernova- 
driven turbulence in the diffuse interstellar medium 
(jMac Low et al.ll200lh . 
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